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ANALYTICAL CORRELATION OF CENTRIFUGAL COMPRESSOR DESIGN 
GEOMETRY FOR MAXIMUM EFFICIENCY WITH SPECIFIC SPEED 

by Michael R. Galvas 
Lewis Research Center 

SUMMARY 

Centrifugal compressor performance was examined analytically to determine op- 
timum geometry for various applications as characterized by specific speed. Seven spe- 
cific losses were calculated for various combinations of inlet tip-exit diameter ratio, in- 
let hub-tip diameter ratio, blade exit backsweep, and inlet tip absolute tangential veloc- 
ity for solid body prewhirl. The losses considered were inlet guide vane loss, blade 
loading loss, skin friction loss, recirculation loss, disk friction loss, vaneless diffuser 
loss, and vaned diffuser loss. Maximum total efficiencies ranged from 0.497 to 0. 868 
for a specific speed range of 0. 257 to 1. 346. 

Curves of rotor exit absolute flow angle, inlet tip -exit diameter ratio, inlet hub -tip 
diameter ratio, head coefficient, and blade exit backsweep are presented over a range of 
specific speeds for various inducer tip speeds to permit rapid selection of optimum com- 
pressor size and shape for a variety of applications. 


INTRODUCTION 

Centrifugal compressors are suitable for aircraft, space-power systems, and other 
applications where small, compact power systems are a necessity. Because the overall 
geometry for each application could vary greatly, it would be advantageous to have a 
method whereby the optimum geometry could be systematically selected. One parameter 
that roughly correlates efficiency with application is specific speed, a term which has 
traditionally been used in the analysis of pumps and compressors. 

Specific speed is a term which relates compressor rotative speed, flow rate, and 
ideal enthalpy rise. Impellers associated with a large total head and low flow rate are 
characterized by a low specific speed and, conversely, those with a low total head and 
high flow rate are characterized by a high specific speed. Many combinations of 



velocity diagram characteristics and compressor geometric parameters are possible at 
a particular specific speed. The combinations of greatest interest are those that will re- 
sult in the maximum attainable efficiency at a particular specific speed. 

Several previous studies have investigated specific speed and its relation to attain- 
able efficiency. Balje (ref. 1) presented plots of efficiency as a function of specific 
speed and specific diameter. Specific diameter is a term relating impeller exit diam- 
eter, adiabatic head, and inlet volume flow rate. From his study it was concluded that, 
for centrifugal compressors with swirl -free inlets, maximum efficiency could be ex- 
pected for backswept impellers with specific speeds in the range of 0. 705 to 1. 018. The 
optimum exit-inlet diameter ratio was mainly a function of specific diameter. Maximum 
attainable efficiency decreased with decreasing specific speed because of gas density 
changes related to high pressure ratios encountered in low specific speed machines. The 
additional losses incurred are due primarily to increased skin friction and disk friction. 

Nichols, McPherson, and Balje (ref. 2) investigated specific -speed - efficiency char- 
acteristics of centrifugal pumps through an analytical determination of losses. Their 
principal findings were that maximum efficiency increased with increasing specific diam- 
eter, that maximum efficiency was not particularly sensitive to impeller inlet hub-tip ra- 
dius ratio, and that maximum efficiency could be expected for backswept impeller blades. 

The general expression of specific speed may be expanded so that several of the 
geometric parameters of interest appear explicitly in the equation. Examination of the 
equation shows that many combinations of the parameters are possible at a given specific 
speed. The combinations of greatest interest are those that will result in maximum at- 
tainable total efficiency at a particular specific speed. The present study is an analyt- 
ical correlation of maximum total efficiency with specific speed by a calculation of spe- 
cific losses that are dependent on geometric parameters and velocity diagram character- 
istics. The method is a one-dimensional analysis conducted at fixed inlet conditions of 
gas state and inducer tip speed. All calculations were made for air entering the com- 
pressor at U.S. standard sea-level conditions. Geometric parameters and velocity char- 
acteristics considered in the study were inlet prewhirl, impeller inlet absolute critical 
velocity ratio, inlet tip -exit diameter ratio, blade backsweep at the impeller exit, and 
inlet hub-tip diameter ratio. Losses considered were those resulting from inlet guide 
vanes if prewhirl was incorporated into the impeller design, impeller blade loading, skin 
friction, recirculation of the working fluid between the vaned diffuser and impeller, dif- 
fusion in the vaneless space, and diffusion in the vaned diffuser. 

Results of this study show the values of geometric variables and velocity diagram 
characteristics that result in maximum efficiency at a given specific speed in the range 
investigated. Individual losses are presented to indicate the changes in loss distribution 
with specific speed along the curves of maximum total efficiency. 
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METHOD OF ANALYSIS 


Approach 

Specific speed, a term relating shaft rotative speed, inlet flow volume, and ideal 
work is defined as 


N 


s 


3/4 


id 


( 1 ) 


Compressor station numbers are shown in figure 1, and all symbols are defined in appen- 
dix A. Equation (1) can be expanded into 



which shows that a number of combinations of geometric variables and velocity diagram 
parameters are possible at a given specific speed. Thus, the problem is to determine 
the combination at a specific speed that will yield maximum efficiency. Individual losses 
were computed for each combination of independent variables considered. The extensive 
range of these combinations results in a variety of different configurations at a particular 
specific speed. Losses considered were those due to inlet guide vanes, blade loading, 
skin friction, disk friction, recirculation, vaneless space diffusion, and vaned diffuser 
diffusion. At the inlet triangles were calculated from a specified absolute critical veloc- 
ity ratio and amount of solid-body prewhirl at the inducer tip (see fig. 2). Gradients in 
axial component of absolute velocity were calculated to satisfy radial equilibrium. At the 
exit relative velocity was assumed to deviate from the blade with a slip velocity calcu- 
lated by the method of reference 3. With the magnitude of relative velocity known from a 
constraint on relative velocity ratio, impeller velocity known from exit -inlet tip diameter 
ratio, specified impeller inlet tip speed, and blade exit angle, the exit triangle was then 
calculated. 

Impeller actual aerodynamic work imparted to the fluid was calculated with the rms 
inlet prewhirl velocity and blade speed and the exit absolute tangential velocity and blade 
speed. Actual shaft work input was determined by adding disk friction and recirculation 
losses to the impeller actual aerodynamic work. The disk friction and recirculation 
work was arbitrarily assumed to escape the system through heat transfer. Isentropic 
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work, which is the work associated with the actual pressure ratio attained, was calcu- 
lated as the actual work minus the summation of internal losses. Adiabatic efficiency 
was calculated as isentropic work divided by actual work. 

Impeller exit geometry was determined from actual aerodynamic work, internal im- 
peller losses, and the constraint on inducer tip-exit relative velocity ratio. 


Calculated Losses 

Inlet guide vane loss. - Equations relating working fluid kinetic energy level, bound- 
ary layer, and blade geometry were developed for axial turbine stators and rotors in ref- 
erence 4. It was assumed that the term for turbine stators was applicable to inlet guide 
vanes of a centrifugal compressor, and the relevant component of the total loss was used 
for the guide vane calculations. This equation was evaluated at mean inlet conditions: 

Ah IGV = 0 ~ 2 ~~ ( V 0 + VlMF ) (3) 

2Re • cot a MF 

Blade loading loss. - Boundary -layer growth in the impeller is highly dependent on 
the internal diffusion of the working fluid. It has been found (ref. 5) that low loss blade 
loadings are those that are characterized by rapid deceleration in the inducer and highly 
loaded inducers in comparison to the outlet. These criteria tend to minimize losses 
associated with secondary flows and separation in the impeller. Reference 6 proposed an 
equation for calculating the diffusion factor of the impeller based on a uniform velocity 
loading along the blade chord. With the diffusion factor calculated by this method, the 
loss due to blade loading was expressed as 

Ah BL = 0. 05d 2 u | (4) 

Skin friction loss. - In addition to the losses resulting from flow separation and sec- 
ondary flows, the impeller incurs losses due to skin friction of the wetted area of the im- 
peller and shroud. Reference 6 proposed an equation for this loss based on pipe friction 
for fully developed turbulent flow: 


Ah gF 


5. 6C f 
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D hyd \ u 2/ £ 
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Disk friction loss . - Windage loss on the compressor back -face was calculated using 
the equation 


Ah DF = °- 01356 


P2 


wRe 


0.2 


4°l 


( 6 ) 


which is a form of the disk friction power loss found in reference 7. 

Recirculation loss . - Losses incurred through the additional work done on the work- 
ing fluid due to back -flow from the vaned diffuser to impeller was expressed as 


Ah RC = u 2 °- 02 (tan «2 )1//2 D f 


(7) 


This is a modification of the loss estimate proposed in reference 6. 

Vaneless diffuser loss . - Between the rotor and vaned diffuser a vaneless space was 
used to decelerate the flow to Mach 0. 8. This was done as a precaution against shock 
flow in the diffuser vanes. Flow angle and Mach number variation with increasing radius 
were determined by numerical solution of the governing differential equations developed 
in reference 8. Total pressure distribution was calculated using the equation developed 
in reference 6. When fluid properties had been determined at the vaned diffuser leading- 
edge radius the vaneless diffuser loss was calculated from the equation: 


Ah VLD ~ C p T 2 


' P 3\ (y l)/7 ( P 3 ' (7 ' 1)/y ' 

vP3/ 


( 8 ) 


Vaned diffuser loss . - Curves of maximum static pressure recovery coefficient at a 
given area ratio were extrapolated from the test data for square throat diffusers reported 
in reference 9 for various combinations of vaned diffuser throat Mach number, blockage, 
and area ratio. Vaned diffuser area ratio and pressure recovery coefficient C** corre- 
sponding to an exit Mach number of 0. 2 were established through iteration. With the exit 
static and total pressures established by this method, vaned diffuser head loss was calcu- 
lated by 


^VD " C p T 2 


'p 4 \ (r - 1)/r fa\b-Wv\ 


V P 4, 




(9) 
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Procedure 


The initial set of performance calculations was made at a fixed inlet absolute critical 
velocity ratio of 0. 3 and inducer tip critical velocity ratio of 0. 520. Subsequent calcula- 
tions were made at inlet absolute critical velocity ratios of 0. 2, 0. 4, and 0. 7 with corre- 
sponding inducer tip critical velocity ratios of 0.346, 0.693, and 1.212, respectively. 
These additional calculations were made to study the effects of internal velocity level on 
compressor performance. For all the calculations, inlet total pressure, total temper- 
ature, total dynamic viscosity, gas constant, and ratio of specific heats were held con- 
stant. All calculations were made for air as the working fluid at U. S. standard sea-level 
conditions. Inducer tip-exit relative velocity ratio was held constant at 0. 7. Low values 
of relative velocity ratio may result in highly separated flow in the rotor, in which case 
the empirical loss estimates would be extremely questionable. Maximum inlet absolute 
tangential velocity for the prewhirl cases was limited to 0.4 of the inducer tip speed. 
Values greater than 0. 4 resulted in inducer velocity diagrams that were not of interest in 
this study. 

The geometric parameters and velocity diagram characteristics d 1 T / D 2> %>’ 
VlT^cr’ anC * V u 1 1 were varied over the following ranges: 
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Inducer tip -exit diameter ratio was varied between 0. 3 and 0. 7 to attain the range of 
specific speeds investigated. Minimum inducer hub-tip diameter ratio was set at 0.3 
since manufacturing and blockage problems are encountered with low values of this var- 
iable. Maximum inducer hub-tip diameter ratio was set at 0. 7. Impeller exit blade back- 
sweep range of 0° to 45° included the maximum efficiency configuration at every specific 
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speed. Vaned diffuser inlet Mach number was limited to 0. 8. A conservative subsonic 
limit was placed on the diffuser inlet Mach number in anticipation of shock-free vaned 
diffuser performance. Vaned diffuser exit Mach number was set at 0. 2 since this repre- 
sents a desirable combustor inlet velocity. 

Since efficiency and specific speed were of primary interest among the output var- 
iables, the calculations were evaluated in these terms. The broad range of input var- 
iables resulted in a wide spread in specific speeds and an appreciable spread in efficien- 
cies at each specific speed. The upper envelope of curves of efficiency as a function of 
specific speed therefore represents the curve of optimum geometry. After this curve had 
been determined, calculated points on or near it were examined for definition of optimum 
geometry. 


RESULTS 

Efficiency Characteristics 

Rotor inlet tip speed was held constant at a critical velocity ratio of 0. 520 for the 
base set of performance calculations with air entering the compressor at U. S. standard 
sea -level conditions. Subsequent calculations were made at rotor inlet tip critical veloc- 
ity ratios of 0.346, 0.693, and 1.212. The combinations of input variables resulted in 
maximum total efficiencies ranging from 0. 497 to 0. 868 and specific speeds from 0. 257 
to 1. 346. Figure 3 shows the variation in maximum total efficiency with specific speed 
and rotor inlet tip speed. The geometric parameters that resulted in these curves were 
examined for definition of optimum geometry. The wide range of specific speeds attained 
with the subsonic inducer tip speeds was not attained with the transonic inducer tip speed 
because the maximum inducer tip-exit diameter ratio for this configuration was limited 
to 0. 5. This results in an impeller exit tip speed of 752 meters per second. Higher 
values were considered to be unreasonable for this study. The large difference in attain- 
able maximum total efficiency for the transonic inducer in comparison with the subsonic 
inducers is due primarily to losses associated with the vaneless and vaned diffusers. 


Compressor Geometry and Velocity Diagram Characteristics 

Compressor geometry was examined for calculated points on or near the curve of 
maximum total efficiency. Several design parameters were then plotted as functions of 
specific speed to reflect the changes in optimum compressor shape with specific speed. 
Figure 4 shows the variation in rotor exit absolute flow angle for a swirl-free inlet 
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Figure 5. - Variation in optimum rotor exit absolute flow angle with specific speed at various inducer tip speeds. 


and four values of solid-body prewhirl for the inducer tip critical velocity ratio of 0. 520. 

The variation in optimum rotor exit absolute flow angle with specific speed and in- 
ducer tip speed is shown in figure 5. The optimum rotor exit absolute flow angle gener- 
ally increased at a given specific speed with increased inducer tip speeds. This resulted 
from higher values of optimum backsweep associated with higher inducer tip speeds. 

Figure 6 shows the variation in optimum inducer hub -tip diameter ratio with specific 
speed for the four inducer tip speeds considered. For the inducer tip critical velocity 
ratio of 0.520, a linear decrease in hub-tip diameter ratio is exhibited in the range of 
specific speeds of 0. 41 to 0. 50. At the inducer tip critical velocity ratio of 0.346 this 
same decrease is noted in the specific speed range of 0. 31 to 0. 39. At the inducer tip 
critical velocity ratios of 0. 693 and 1. 212 the decrease is exhibited in the specific speed 
ranges of 0. 59 to 0. 68 and 0. 79 to 1. 01, respectively. 



Figure 6. - Variation in optimum inducer hub-tip diameter ratio with specific speed at various inducertip speeds. 
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The optimum inlet tip-exit diameter ratio is shown in figure 7. The low values at 
low specific speeds reflect the high ideal work relative to inlet flow volume that produces 
low specific speeds. This ratio increases linearly to the specified maximum at a specific 
speed of 0. 95 for the inducer tip critical velocity ratios of 0. 36 and 0. 520. This char- 
acteristic curve was translated to higher specific speed values for the inducer tip critical 
velocity ratios of 0. 693 and 1. 212. 
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The variation in blade exit backsweep with specific speed for the four inducer tip 
speeds is shown in figure 8. Backsweep has a profound influence on impeller exit abso- 
lute Mach number reduction for a fixed -reaction rotor, so larger values of backsweep 
would be expected at lower specific speeds where higher absolute Mach numbers increase 
vaneless diffuser losses. However, for the impeller tip critical velocity ratio of 0. 346 
exit Mach numbers were near 0. 8 maximum. This resulted in good vaned diffuser per- 
formance, and the improvements in diffuser pressure recovery attained with backsweep 
were negated by increased recirculation losses in the lower specific speed region where 
rotor exit absolute flow angles are large. In the higher specific speed range moderate 
backsweep was beneficial because recirculation losses were not significantly increased. 
For the inducer tip critical vebcity ratio of 0. 520 backsweep was beneficial up to spe- 
cific speeds of 0. 59. For the two higher inducer tip critical velocity ratios the impeller 
exit absolute Mach number was always greater than 0. 8 and backswept impellers pro- 
duced maximum total efficiency over the entire range of specific speeds investigated. 
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Figure 10. - Variation in optimum rotor exit blade height -diameter ratio with specific speed at various inducer tip 
speeds. 


In the absence of a transonic loss correlation, compressor configurations with swirl - 
free inlets resulted in maximum total efficiency over the range of specific speeds investi- 
gated. If a correlation reflecting shock - boundary -layer interaction were included for 
configurations with transonic inducers, prewhirl would be expected to increase efficiency. 
However, prewhirl reduces the attainable pressure ratio for a given rotative speed and 
inlet tip-exit diameter ratio. 

The ratio of exit blade height -diameter for maximum total efficiency is shown in fig- 
ure 9 for a swirl -free inlet and four values of solid -body prewhirl at the inducer tip crit- 
ical velocity ratio of 0. 520. Increasing ratios with increasing specific speed reflect gas 
density decreases with compressor pressure ratio decreases. The effect of inducer tip 
speed on exit blade height -diameter ratio for the four configurations with swirl-free inlet 
only is shown in figure 10. Blade height decreases with increasing tip speed. This de- 
crease is due to higher gas density and compressor pressure ratio associated with in- 
creased rotative speed. 

Figure 11 shows the variation in compressor total pressure ratio for the four inducer 
tip speeds considered. The high-speed configuration had a transonic relative gas velocity 
at the inducer tip, and in the absence of transonic rotor losses the predicted attainable 
pressure ratio is obviously optimistic. 
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Figure 11. - Variation in compressor total pressure ratio for maximum total efficiency configuration with specific speed 
at various inducer tip speeds. 


Loss Distribution 

The distribution of internal losses along the curves of maximum total efficiency as a 
function of specific speed is shown in figure 12 for the four inducer tip speeds studied. 
The changes in loss distribution are shown to indicate the effect of internal velocity level 
on compressor performance. In the low range of specific speed attained by each inducer 
tip speed, the losses are dominated by disk friction and vaneless and vaned diffuser 
losses. Disk friction losses are large in the lower specific speed range because of the 
higher exit tip speeds and gas density that characterize low specific speed machines. 
Vaneless and vaned diffuser losses are large because of the higher impeller exit absolute 
Mach numbers associated with the lower specific speed machines. Skin friction loss, 
which was dependent on inducer inlet hub -tip diameter ratio, showed a marked change in 
level in the specific speed range where optimum inducer hub-tip diameter ratio changed 
for each tip speed. 
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Figure 12. - Variation in internal loss distribution along curve of maximum total efficiency with specific speed. 








Figure 12. - Concluded. 


Head Coefficient 


The variation in head coefficient ^H^lu^ with specific speed and inducer tip speed 
is shown in figure 13. The variation in this parameter is included in addition to figure 11 
because it permits a straightforward determination of impeller exit tip speed for a de- 
sired pressure ratio. The decrease in head coefficient in the mid-specific -speed range 
for the inducer tip critical velocity ratio of 0. 346 is associated with the rapid change 
from radial to backswept impellers for maximum total efficiency. Efficiency is insen- 
sitive to backsweep in comparison to actual impeller work, thereby reducing head coeffi- 
cient which is dependent on isentropic work. 



Specific speed, 


Figure 13. - Variation in optimum head coefficient with specific speed at various inducer tip speeds. 


Application of Results 

The results of this investigation provide a basis on which centrifugal compressors 
can be sized for applications in which maximum total efficiency is an important design 
parameter. Compressor specific speed can be estimated from desired rotative speed, 
pressure ratio, weight flow, and inlet density. When the approximate specific speed is 
known, figure 13 can be used to determine the rotor exit speed and, consequently, rotor 
exit diameter, and figures 6 and 7 can then be used to size the inlet. Then a new estimate 
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of specific speed can be made from the indicated volume flow rate required for the weight 
flow and inlet dimensions indicated. Several iterations of this procedure will determine 
compressor size for the application considered. Figures 4, 8, and 9 can then be used to 
determine impeller exit velocity diagram characteristics and geometry. Caution should 
be used with figures 6 and 9 because mechanical considerations are not reflected in these 
graphs. Severe inlet blockage problems may be encountered with low inlet hub-tip diam- 
eter ratios indicated in the high specific speed range. Actual exit blade height should be 
calculated from continuity at the rotor exit. The rotor exit velocity diagram can be con- 
structed with the aid of figure 4, the wheel speed obtained from the use of figure 13, and 
the constraint placed on impeller exit -inlet tip relative velocity ratio. Impeller exit total 
temperature can be calculated from the actual aerodynamic work indicated by the impel- 
ler exit velocity triangle. Impeller exit total pressure can be estimated using a rotor 
efficiency computed with the aid of figure 12. Determination of the impeller exit absolute 
Mach number permits calculation of the exit static density. True blade height can then 
be calculated from the desired weight flow, the calculated static density, the indicated 
impeller exit diameter, and the absolute meridional velocity. 

Hub and shroud contours and blade shape resulting in desirable diffusion may then be 
established by use of an internal velocity calculation program such as reference 10. 


SUMMARY OF RESULTS 

Centrifugal compressor geometry was examined analytically to determine the values 
of geometric variables that yield maximum total efficiency over a range of specific 
speeds. Geometric parameters considered were inlet hub-tip diameter ratio, inlet tip- 
exit diameter ratio, rotor exit blade height-exit diameter ratio, and blade exit back- 
sweep. The variation of these values with specific speed was then determined, thereby 
providing a method by which low -loss geometry can be rapidly selected for any applica- 
tion in the range of specific speed investigation. The results of the analysis can be sum- 
marized as follows: 

1. The flow conditions and geometric variables resulted in a specific speed range of 
0. 257 to 1. 346 for the inducer tip speeds studied. Maximum total efficiencies ranged 
from 0.497 to 0. 868. 

2. Efficiency dropped rapidly with decreasing specific speeds in the lower specific 
speed range attainable by each inducer tip speed. This drop resulted from large losses 
associated with disk friction and both vaneless and vaned diffusers. 

3. All of the geometric variables and velocity diagram characteristics studied were 
sensitive to changes in inducer tip speed. 
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4. Swirl-free inlet configurations resulted in maximum total efficiency in the range 
of specific speeds investigated. This resulted from a neglect of shock losses in the 
impeller. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 21, 1972, 

132-15. 


19 



APPENDIX A 


SYMBOLS 


a local acoustic velocity, m/sec 

B diffuser effective depth ratio 

b blade height, m 

Cj skin friction coefficient 

C specific heat, J/(kg)(K) 

r 

C* * maximum pressure recovery coefficient at a given area ratio 

r 

D diameter, m 

Dj diffusion factor 

e g inlet guide vane loss coefficient 

H overall ideal compressor enthalpy, J/kg 

Ah incremental compressor enthalpy, J/kg 

L blade length, m 

M Mach number 

N g specific speed, dimensionless 

2 

p pressure, N/m 

Q volume flow rate, m /sec 
q dimensionless enthalpy 

R gas constant, J/(kg)(K) 

R radius ratio 

Re Reynolds number 

r radius , m 

T temperature 

u blade speed, m/sec 

V absolute gas velocity, m/sec 

W relative gas velocity, m/sec 

w mass flow rate, kg/sec 
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number of blades 

absolute flow angle, deg from meridional 
relative flow angle, deg from meridional 
blade angle, deg from meridional 
ratio of specific heats 
deviation angle between flow and blade 
incremental boundary -layer displacement thickness, 
loss coefficient 
efficiency 

decrement in efficiency 
hub -tip diameter ratio 

2 

dynamic viscosity, N-sec/m 

2 

kinematic viscosity, m /sec 
3 

gas density, kg/m 
flow coefficient 
head coefficient 
impeller angular velocity 
Subscripts: 
act actual 

aero aerodynamic 

av average 

BL blade loading 

calc calculated 

cr critical condition 

DF disk friction 

est estimated 

H hub 

hyd hydraulic 

IGV inlet guide vane 

i iteration subscript 


Z 

a 

J3 

% 

r 

5 

A 5 

5 

V 

Atj 

X 

M 

v 

P 

<P 

(x> 



id ideal 

MF root-mean-square 
m meridional 

R rotor 

RC recirculation 
SF skin friction 

SL slip 

T tip 

th theoretical 

u tangential 

VD vaned diffuser 

VLD vaneless diffuser 

0 station just upstream of inlet guide vanes 

1 impeller inlet 

2 impeller exit 

3 vaned diffuser inlet 

4 vaned diffuser exit 

Superscripts : 

absolute stagnation 
relative stagnation 
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APPENDIX B 


EQUATIONS 

The equations used in the calculation of compressor performance are listed here in 
the general order of solution. Several quantities were held constant for each set of cal- 
culations. These constants were inlet total conditions of the working fluid Pq and Tq, 
fluid properties y, R, and and inducer tip speed ju^p. 

Geometric parameters D^/Dg, X, Vuj T , and were varied to obtain the range 

of specific speeds. Initial calculations were made at U iq>/V cr °* 52, and further in- 

vestigations were performed at values of 0. 346, 0.693, and 1. 212 to study the effects of 
compressibility on compressor geometry and performance. 


Impeller Inlet Velocity Triangles 

For given values of V^/V^ and prewhirl velocity at the impeller tip, relative 
velocity and meridional component of absolute velocity were calculated. With the solid 
body vortex prewhirl considered in this analysis, tangential components of absolute 
velocity were calculated by the relations 


Vu 


IMF 




(Bl) 


Vu 


1H 



(B2) 


where D^p is the root -mean-square inlet diameter. Meridional component of abso- 
lute velocity must vary according to the relation (ref. 11) 


Vmrf 7 ^ < B3) 

in order to satisfy radial equilibrium in a solid-body vortex. The constant in equa- 
tion (B3) was evaluated at the tip since Vm^ had been determined by the specified in- 
let critical velocity ratio and tip prewhirl velocity. With the constant known, absolute 
meridional components of velocity were calculated at the mean and hub. The inlet tri- 
angles were then determined from known wheel speeds and absolute velocity components. 
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Impeller Exit Velocity Diagram 


At the impeller exit the flow was assumed to deviate from the blade with a slip veloc- 
ity calculated by the equation from reference 3 


SL 


^COS P 


2b 


Ur 


-0.7 


(B4) 


where Z is the number of blades proposed by reference 12 


v R R °2 + D 1 lav + ^2b’ 

Z - 6. 5 cos' 

D 2 " D 1 

where 


(B5) 


D, = ^0. + Dj h ) 

The deviation angle between the blade and the flow is then specified by 


(B6) 


S = sin-1 V S L ^ 

w„ 


(B7) 


The other velocity diagram characteristics (see fig. 2) were calculated from the 
relations 


Vm 2 = W 2 cos (/3 2b + 6) 

V u2 = u 2 " V SL " W 2 cos ^2b + 6 ^ tan ^2b 


(B8) 

(B9) 


V ' 


°2 


= tan 


-1 


+ Vm 2 

(BIO) 

/Vu„\ 



(Bll) 

V m 2/ 
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u 2 ~ U 1T 


Inlet Guide Vane Loss 


This loss was calculated using the equation 


2Re 4 * cot a* 


y2 + 

V 0 IMF, 


(B13) 


Total pressure loss across the guide vanes was calculated by the method of reference 13, 
assuming that turbine stator loss coefficients could be applied to centrifugal compressor 
inlet guide vanes. The guide vane loss coefficient was calculated from the equation 


0. 0076 

cos a -| 1 ^ " 0. 025 


Ideal kinetic energy leaving the vanes was calculated by 


(B14) 


<KE) irt = 


(B15) 


where 


KE = 


(B16) 


Mean static temperature downstream of the inlet guide vanes was calculated from 
the equation 


T 1MF ~ T 0 


(B17) 
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The remaining fluid state conditions were computed from the relations: 



y/(v- 1) 


IMF! 


2C P T 0, 


P 1MF 

p i 


(KE) id\ r/<y " 1) 

vy 


p 


IMF 



p 0 


p 


IMF " P 1MF 



IMF 


P 1MF 


? 

P 1MF 


t 



’ (A 1/y 

p IMF " p lMF “ T 

\ p /IMF 


P 0 


RT' 


Velocity upstream of the inlet guide vanes was calculated by trial and error 
of the equation 


p 0 V 0 

p 0 V cr 


1 - 


r- i ( v o 

y + 1 \V 


cr^ 


i/(y-D 


cr 


(B18) 


(B19) 


(B20) 


(B21) 


(B22) 


(B23) 


(B24) 

solution 


(B25) 
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where 


P0 V 0 ~ P 1MF V 1MF COS a lMF 


(B26) 


Blade Loading Loss 


The work loss due to blade loading was previously expressed as 


ih BL = o. 05 d|u| 


(B27) 


as suggested by reference 6. Where the diffusion factor Dj was defined as 


0. 75 


D f = 0. 3 + 


% 


W 


IT 


Ur 


5(l-^ + 2 DlT 


7T 


Do 


D r 


(B28) 


for W 2 /W 1T =0.7 

The dimensionless actual head was obtained from blade exit speed and absolute tan- 
gential velocity and rms inlet wheel speed and prewhirl velocity: 


**th 


U 2 Vu 2 ‘ u lMF Vu lMF 


Ur 


Skin Friction Loss 


Skin friction loss was found from the equation of reference 6: 


Ah - 5 6C M 2 u 2 

Ah gF -5.6C f7 - v (— ] 2 


D hyd\ \ u 2 / 


av 


Dr 


(B29) 


(B30) 


where the mean flow path blade length -to -exit diameter ratio was expressed as 
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1 


IMF 


L 1 0. 3048 


Dg 2 cos 0 


2b 


(B31) 


The ratio of mean hydraulic -to -exit diameter was found from the equation 
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IT 


°hyd _ 


D n 


D 0 


D r 


2Z 


u cos 0 


2b u 2 1 - X jt(1 + X) 


V 


1 + ten2 *1T 


(B32) 


and the ratio of mean relative velocity -to -impeller exit velocity was determined by 


r 


/w 

\u„ 




'Vm 


IMF! 


av 


u. 




since W 2 /Wj T = 0. 7 


+ +0.49 

Do, 


Vm 


IMF 


"2 


2 \2 

. IT 

+ 

D„ 


(B33) 


Recirculation Loss 


The recirculation loss was simply a function of the exit absolute gas angle and im- 
peller diffusion factor (ref. 6): 


Ah RC = 0. 02(tan ct 2 )^ 2 D^u 


(B34) 


Disk Friction Loss 

Impeller disk friction loss was computed with the equation of reference 7: 


Ah DF = 0.01356 


P 2 


wRe 


0.2 


u 2 D 2 


(B35) 
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in which the weight flow rate w was calculated from the mean flow properties and inlet 
geometry 


w “ “ D 1T ^ " a2 ^1MF V 1MF cos 01 IMF ( B36 ^ 

Reynolds number for this loss calculation was based on impeller exit dimensions and in- 
let stagnation conditions : 


He = U2l>2PMF (B37: 

Mo 

Calculation of the exit density required iterative computations. Exit density corre- 
sponding to an isentropic pressure rise across the rotor was chosen as a first approx - 
imitiation and the rotor adiabatic efficiency, which reflects blade loading, skin friction, 
and inlet guide vane losses, was calculated using the estimated value. The equations 
used in this procedure were as follows: 
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T 1MF - T 1MF + 


IMF 


2C 


2 2 

_ -r ,T . °2 _U 1MF 

i 2 - 1 1 + 


2C 


P 


1 1 vv o 

T 2 = T 2 - — 
2C 


T 2 T 2 + 


2C 


Ah , = C T n — - l 1 
act p 01 m t 

V X 0 


Ah^ = Ah„_j. - Ah Tmr — Ahn, - Ah c 


act ' IGV “ BL ' au SF 


(B38) 

(B39) 

(B40) 

(B41) 

(B42) 

(B43) 
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(B44) 



Impeller Exit Geometry 

Impeller exit blade height was calculated with the aid of the continuity equation. 

Blade height necessary to pass the calculated mass flow is 

b 9 = ^ (B48) 

71 P 2°2 V m2 

Vaneless Diffuser Loss 

Vaneless diffuser radius ratio was determined by the restrictions on the Mach num- 
ber entering the vaned diffuser. Vaned diffuser inlet Mach number was limited to 0. 8 to 
guard against shock formation in the vaned diffuser. Because of the wide range of var- 
iables considered, the rotor exit absolute Mach number varied extensively. In cases 
where its value was less than or equal to 0. 8, the vaneless space radius ratio was arbi- 
trarily set at 1.02; in cases where its value was greater than 0. 8, the radius ratio was 
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determined by the amount of diffusion necessary to reduce the Mach number to the spec- 
ified limit. 

Reference 8 developed differential equations relating Mach number, flow angle, and 
total temperature to vaneless diffuser radius ratio through the fundamental relations of 
continuity, meridional equilibrium, tangential equilibrium, heat transfer, and fluid state. 
Meridional equilibrium was satisfied by a balance of pressure forces, shear forces, and 
inertia forces. Tangential equilibrium was satisfied by a balance of shear forces and 
inertia forces. Combination of the equilibrium equations with continuity and state equa- 
tions resulted in the equations relating the changes in Mach number and flow angle to 
changes in radius ratio. 

Mach number and flow angle were calculated at increments of 0. 02 in radius ratio by 
the numerical method described in detail in reference 8. The governing differential 
equations for adiabatic diffusion in a radial passage are 


and 


1 dM" 


- 2^1 + 


^1m 2 
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2 dR 


M 2 - sec 2 a 


(■} M 2 - tan 2 a) s 


1 DB sec' 1 a 


B cos a B dR 
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(B49) 


1 d(tan a) _ sec z a 


tan a dR 


M 2 - sec 2 a 


1 + (y - 1)M" 


B cos a 


1 dB_ M 
B dR R 


(B50) 


where 


s = 


C f r 2 


R = . 


B = — 
b 0 


(B51) 

(B52) 

(B53) 


The schedule of effective passage depth was determined by boundary -layer displace- 
ment thickness growth on the end walls. The flow was assumed to follow a log-spiral 
path determined by the calculated angle along each increment in radius ratio. The incre- 
mental displacement thickness at each radius ratio was then calculated from the equation 
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-1/5 


(B 54) 


* 

A6 


0.37 ARr„ /V 0 ARr, 


8 cos a \ v cos a , 


which is the displacement thickness predicted by the one-seventh power velocity distribu- 
tion law discussed in reference 14. 

The effective passage depth was then 


B i + 1 = B i 


2 A6 


(B55) 


and 


AB = B i - B i + 1 (B56) 

Loss in total pressure in the vaneless space was computed from the equation found 
in reference 5: 

3 ©(^) 5dS 

(B57) 

The integral in the above equation can readily be determined by numerical methods 
using the Mach number distribution calculated by equation (B49) and the familiar isen- 
tropic flow relations of reference 15. The radius ratio at which the vaned diffuser inlet 
Mach number restriction was satisfied was also assumed to be vaned diffuser throat. 
Throat static pressure was calculated from the equation 


= 1 + 


yC 


f r 2 
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cos at 2 b 2 M 
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Vaneless diffuser loss was then determined by 


2 \y/ (y- 1) 

2 v 


Ah VLD ~ C p T 2 


'P 8 \ (y ' 1)/y /P3 x(r_1)/y 


, p 3, 


V P 2y 
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(B59) 
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Vaned diffuser performance was predicted from test data reported in reference 9. 
Lines of maximum pressure recovery for a prescribed area ratio were estimated from 
the performance maps reported for single plane divergence diffusers with square throats. 
Throat conditions of total pressure, Mach number, and blockage were determined by the 
analysis of the vaneless diffuser. Vaned diffuser area ratio required to decelerate the 
flow to an exit Mach number of 0. 2 was determined by iteration. Using the one- 
dimensional isentropic area ratio 


M, 

AR = — 

0.2 


1 + 21 


1 + 


— ( 0 . 2 ) 

2 


2H-I M? 


2 


-,(y+l)/2(y-l) 


(B60) 


as an initial guess, the maximum pressure recovery coefficient Cp corresponding to 
the throat Mach number, blockage, and area ratio was interpolated. Exit static and total 
pressures were then calculated from the pressure recovery coefficient and isentropic 
relation, respectively: 


^ ^ f 

P 4 = c p (P 3 -P 3 ) + P 3 


(B61) 


, r v 1 oiy/(y-D 
p 4 = P 4 L 1 + (o a J 


(B62) 


The subsequent approximations to the actual required area ratio are then a function 
of the previous approximation and the estimated and calculated exit total pressures: 


AR. 


p 4 i 

i+1 = AR i 

P4.i+1 


(B63) 


The iteration was continued until two consecutive estimates of area ratio agreed 
within 0.001. When convergence was achieved, vaned diffuser loss was determined from 
the equation: 


Ah VD ~ C p T 2 


'P4 


(y-l)/y 



(B64) 
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Overall Efficiency 


The overall adiabatic efficiency was calculated as follows: 


n . „ 


Ah a< 


- f Ah 


(B65) 


Efficiency Decrements 

Decrements in efficiency caused by the individual losses were obtained as follows : 
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where 


h act ^aero + 


Ah 


DF 


Ah 


RC 


Specific Speed 


This term is defined as 


and is computed using 




N s = 


V^i 


H. 


3/4 

id 


2u 


u> = 


IT 


D 


IT 


H id _ t? Ad( u 2 Vu 2 “ u 1MF Vu 1Mf) 
Ql = ~ D 1T ^ " x2 ) V 1MF cos “IMF 


(B73) 


(B74) 

(B75) 

(B76) 

(B77) 


35 



REFERENCES 


1. Balje, O. E. : A Study on Design Criteria and Matching of Turbomachines: 

Part B - Compressor and Pump Performance and Matching of Turbo components. 

J. Eng. Power, vol. 84, no. 1, Jan. 1962, pp. 103-114. 

2. Nichols, K. E.; McPherson, D. G. ; and Balje, O. E. : Study of Turbine and Turbo - 

pump Design Parameters. Rep. S/TD 1735, Turbo Div, Sundstrand Corp. , 

Nov. 20, 1959. 

3. Wiesner, F. J. : A Review of Slip Factors for Centrifugal Impellers. J. Eng. 

Power, vol. 89, no. 4, Oct. 1967, pp. 558-572. 

4. Stewart, Warner L. : A Study of Axial-Flow Turbine Efficiency Characteristics in 

Terms of Velocity Diagram Parameters. Paper 61-WA-37, ASME, Nov. 1961. 

5. Dallenbach, F. : The Aerodynamic Design and Performance of Centrifugal and 

Mixed-Flow Compressors. Technical Progress Series. Vol. 3. SAE, 1961. 

6. Coppage, J. E. ; Dallenbach, F. ; Eichenberger, J. P. ; Hlavaka, G. E. ; 

Knoernschild, E. M. ; and Vanke, N. : Study of Supersonic Radial Compressors 
for Refrigeration and Pressurization Systems. AiResearch Mfg. Co. (WADC TR 
55-257), Dec. 1956. 

7. Shepherd, Dennis G. : Principles of Turbo machinery. Macmillan Co. , 1956. 

8. Stanitz, John D. : One-Dimensional Compressible Flow in Vaneless Diffusers of 

Radial - and Mixed -Flow Centrifugal Compressors, Including Effects of Friction, 
Heat Transfer, and Area Change. NACA TN 2610, 1952. 

9. Runstadler, Peter W. , Jr.: Pressure Recovery Performance of Straight -Channel, 

Single-Plane Divergence Diffusers at High Mach Numbers. Rep. N-88, Creare, 
Inc. (USAAVLABS -TR-69 -56 , AD-865300), Oct. 1969. 

10. Katsanis, Theodore: Use of Arbitrary Quasi -Qrthogonals for Calculating Flow 

Distribution in the Meridional Plane of a Turbomachine. NASA TN D-2 546, 1964. 

11. Vavra, Michael H. : Aero -Thermodynamics and Flow in Turbomachines. 

John Wiley & Sons, Inc. , 1960. 

12. Lazarkiewicz, Stephen; and Troskolanski, A. T. : Impeller Pumps. 

Pergamon Press, 1965. 

13. Rohlik, Harold E. : Analytical Determination of Radial Inflow Turbine Design 

Geometry for Maximum Efficiency. NASA TN D-4384, 1968. 


36 



14. Schlichting, Hermann (J. Kestin, trans.): Boundary Layer Theory. Sixth ed. , 

McGraw-Hill Book Co. , Inc. , 1968. 

15. Shapiro, Ascher H. : The Dynamics and Thermodynamics of Compressible Fluid 

Flow. Vol. 1, Parts 1 and 2. Ronald Press Co. , 1953. 


NASA -Langley, 1972 1 E-6638 


37 




